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Abstract 

A search for events with a high energy isolated electron or muon and missing transverse 
momentum has been performed at the electron-proton collider HERA using an integrated 
luminosity of 13.6 pb _1 in e~p scattering and 104.7 pb _1 in e + p scattering. Within the 
Standard Model such events are expected to be mainly due to W boson production with 
subsequent leptonic decay In e~p interactions one event is observed in the electron chan- 
nel and none in the muon channel, consistent with the expectation of the Standard Model. 
In the e + p data a total of 18 events are seen in the electron and muon channels compared 
to an expectation of 12.4 ± 1.7 dominated by W production (9.4 ± 1.6). Whilst the overall 
observed number of events is broadly in agreement with the number predicted by the Stan- 
dard Model, there is an excess of events with transverse momentum of the hadronic system 
greater than 25 GeV with 10 events found compared to 2.9 ± 0.5 expected. The results 
are used to determine the cross section for events with an isolated electron or muon and 
missing transverse momentum. 
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1 Introduction 



The HERA collaborations HI and ZEUS have previously reported [1-3] the observation of 
events with an isolated high energy lepton and missing transverse momentum in e + p collisions 
recorded during the period 1994-1997. The dominant Standard Model (SM) contribution to this 
topology is real W boson production with subsequent leptonic decay. Such events can also be a 
signature of new phenomena beyond the Standard Model [4]. HI has reported [2] one e~ event 
and 5 ^ events compared to expectations from the Standard Model of 2.4 ± 0.5 and 0.8 ± 0.2 
for the e ± and ^ channels respectively, with W contributions of 1 .65 ±0.47 (e) and 0.53 ± 0. 1 1 
(//). For the same data taking period ZEUS has reported [3] 3 (0) (//^ events compared to 
an expectation of 2.1 (0.8) W events and 1.1 ± 0.3 (0.7 ± 0.2) events from other processes. In 
the present paper a search for events with isolated electrons 1 or muons and missing transverse 
momentum is performed in an extended phase space and with improved background rejection. 
The complete HERA I data sample (1994-2000) is analysed here. This corresponds to an 
integrated luminosity of 118.4 pb _1 , which represents a factor of three increase with respect to 
the previous published result. 

This paper is organised as follows. Section 2 describes the SM processes that contribute to 
the signal and to the background. Section 3 describes the HI detector and experimental con- 
ditions. Section 4 outlines the lepton identification criteria and the reconstruction methods for 
the hadronic final state. The selection requirements for the electron and muon channels are 
described in section 5. Studies of background processes are presented in section 6. Section 
7 deals with systematic uncertainties and section 8 presents the results of the analysis includ- 
ing the numbers of events seen, the kinematics of the selected events and the measured cross 
sections. The results of a search for W production in the hadronic decay channel are given in 
section 9. The paper is briefly summarised in section 10. 

2 Standard Model Processes 

The processes within the Standard Model that are expected to lead to a final state containing 
an isolated electron or muon and missing transverse momentum, due to penetrating particles 
escaping detection in the apparatus, are described in detail in [2] and are only briefly outlined 
in this section. The processes are called "signal" if they produce events which contain a gen- 
uine isolated electron or muon and genuine missing transverse momentum in the final state. 
The processes are defined as "background" if they contribute to the selected sample through 
misidentification or mismeasurement. For the background processes, a fake lepton, fake miss- 
ing transverse momentum or both can be reconstructed and may lead to the topology of interest. 
The following processes are considered. 

W production: ep — > eW ± X or ep — > vW ± X (signal) 

Real W production in electron proton collisions with subsequent leptonic decay W — > Iv, 
proceeding via photoproduction, is the dominant SM process that produces events with 

'in this paper "electron" refers generically to both electrons and positrons. Where distinction is required the 
terms e~ and e + are used. 
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prominent high P T isolated leptons and missing transverse momentum. W bosons are 
predicted to be produced mainly in resolved photon interactions, in which the W typically 
has small transverse momentum, whilst in direct photon interactions the W transverse 
momentum may be larger. 

In this paper, the SM prediction for W production via ep — > eW ± X is calculated by 
using a next to leading order (NLO) Quantum Chromodynamics (QCD) calculation [5] in 
the framework of the EPVEC [6] event generator. Each event generated by EPVEC ac- 
cording to its default LO cross section is weighted by a factor dependent on the transverse 
momentum and rapidity of the W [7], such that the resulting cross section corresponds to 
the NLO calculation. The ACFGP [8] parameterisation is used for the photon structure 
and the CTEQ4M [9] parton distribution functions are used for the proton. The renor- 
malisation scale is taken to be equal to the factorisation scale and is fixed to the W mass. 
Final state parton showers are simulated using the PYTHIA framework [10]. 

The NLO corrections are found to be of the order of 30% at low W transverse momentum 
(resolved photon interactions) and typically 10% at high W transverse momentum (direct 
photon interactions) [5]. The NLO calculation reduces the theory error to 15% (from 30% 
at leading order). 

The charged current process ep — > uW ± X is calculated with EPVEC [6] and found to 
contribute less than 7% of the predicted signal cross section. 

The total predicted W production cross section amounts to 1.1 pb for an electron-proton 
centre of mass energy of y/s = 300 GeV and 1.3 pb for y^i = 318 GeV. 

Z production : ep — > eZ(-^ vv)X (signal) 

A small number of signal events may be produced by Z production with subsequent 
decay to neutrinos. The outgoing electron from this reaction can scatter into the detector 
yielding the isolated lepton in the event while genuine missing transverse momentum is 
produced by the neutrinos. This process is calculated with the EPVEC generator and 
found to contribute less than 3% of the predicted signal cross section. 

Charged Current (CC) processes : ep — > vX (background) 

A CC deep inelastic event can mimic the selected topology if a particle in the hadronic fi- 
nal state or a radiated photon is interpreted as an isolated lepton. The generator DJANGO 
[1 1] is used to calculate this contribution to the background. 

Neutral Current (NC) processes : ep — > eX (background) 

The scattered electron in a NC deep inelastic event yields an isolated high energy lepton, 
but measured missing transverse momentum can only be produced by fluctuations in the 
detector response or by undetected particles due to limited geometrical acceptance. The 
generator RAPGAP [12] is used to calculate this contribution to the background. 

Photoproduction of jets : ^p — > X (background) 

The generator PYTHIA [13] is used to calculate the contribution from hard scattering 
photoproduction processes. Background from this process may occur if a particle from 
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the hadronic final state is interpreted as an isolated lepton and missing transverse mo- 
mentum is measured due to fluctuations in the detector response or limited geometrical 
acceptance. 

Lepton pair (LP) production : ep — ► e l + l~X (background) 

Lepton pair production can mimic the selected topology if one lepton escapes detec- 
tion and measurement errors cause apparent missing momentum. The generator GRAPE 
1.1 [14], based on a full calculation of electroweak diagrams, is used. The dominant con- 
tribution is due to photon-photon processes and is cross-checked with the LPAIR [15] 
generator. Internal photon conversions are also calculated. Z production and its subse- 
quent decay into charged leptons is also included in GRAPE. This contribution is found 
to be negligible. 

In order to determine signal acceptances and background contributions, the detector re- 
sponse to events produced by the above programs is simulated in detail using a program based 
on GEANT [16]. The simulated events are then subjected to the same reconstruction and anal- 
ysis chain as the data. 

3 Experimental Conditions 

Results are presented for the 37.0 pb 1 of e + p data taken in 1994-1997 at an electron-proton 
centre of mass energy of y/s = 300 GeV , the 13.6 pb" 1 of e~p data (1998-1999, = 318 
GeV) and the 67.7 pb" 1 of e + p data (1999-2000, ^1 = 318 GeV). 

A detailed description of the HI detector can be found in [17]. Only those components of 
particular importance to this analysis are described here. 

The inner tracking system consisting of central and forward 2 tracking detectors (drift cham- 
bers) is used to measure charged particle trajectories and to determine the interaction vertex. A 
solenoidal magnetic field allows the measurement of the particle transverse momenta. 

Electromagnetic and hadronic final state particles are absorbed in a highly segmented Liquid 
Argon (LAr) calorimeter [18]. The calorimeter is 5 to 8 interaction lengths deep depending on 
the polar angle of the particle. A lead-fibre calorimeter (SpaCal) is used to detect backward 
going electrons and hadrons. 

The LAr calorimeter is surrounded by a superconducting coil with an iron return yoke in- 
strumented with streamer tubes. Tracks of muons, which penetrate beyond the calorimeter, are 
reconstructed from their hit pattern in the streamer tubes. The instrumented iron is also used as 
a backing calorimeter to measure the energy of hadrons that are not fully absorbed in the LAr 
calorimeter. 

In the forward region of the detector a set of drift chamber layers (the forward muon system) 
detects muons and, together with an iron toroidal magnet, allows a momentum measurement. 
Around the beam pipe, the plug calorimeter measures hadronic activity at low polar angles. 

2 The forward direction and the positive z-axis are taken to be that of the proton beam direction. All polar 
angles are defined with respect to the positive z-axis. 
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The LAr calorimeter provides the main trigger for events with high transverse momentum. 
The trigger efficiency is ~ 98% for events with an electron which has transverse momentum 
above 10 GeV. For events with high missing transverse momentum, determined from an imbal- 
ance in transverse momentum measured in the calorimeter P|. al °, the trigger efficiency is ~ 98% 
when P| al ° > 25 GeV and is ~ 50% when P£ al ° = 12 GeV [19]. Events may also be triggered 
by a pattern consistent with a minimum ionising particle in the muon system in coincidence 
with tracks in the tracking detectors. 



4 Lepton Identification and Hadronic Reconstruction 

An electron candidate is defined [20] by the presence of a compact and isolated electromagnetic 
cluster of energy in the LAr calorimeter, with the requirement of an associated track having an 
extrapolated distance of closest approach to the cluster of less than 12 cm. Electrons found 
in regions between calorimeter modules containing large amounts of inactive material are ex- 
cluded [19]. The energy of the electron candidate is measured from the calorimeter cluster. 
The additional energy allowed within a cone of radius 1 in pseudorapidity-azimuth (rj-cf)) space 
around the electron candidate is required to be less than 3% of the energy attributed to the elec- 
tron candidate. The efficiency of electron identification is established using NC events and is 
greater than 98% [19]. 

A muon candidate is identified by a track in the forward muon system or a charged track 
in the inner tracking system associated with a track segment or an energy deposit in the instru- 
mented iron. The muon momentum is measured from the track curvature in the solenoidal or 
toroidal magnetic field. A muon candidate may have no more than 8 GeV deposited in the LAr 
calorimeter in a cone of radius 0.5 in (rj-cf)) space associated with its track. The efficiency to 
identify muons is established using elastic LP events [21] and is greater than 90%. 

Identified leptons are characterised by the following variables, where I represents e or /i: 

• Pip, the transverse momentum of an identified muon or electron; 

• 6i, the polar angle of the muon or electron. 

In order to check that the probability to misidentify a particle as an electron or muon is 
well described by the simulation, a sample of NC events is used, in which a second electron 
or a muon is found in the event. In the majority of cases this second lepton results from the 
misidentification of a hadron from the final state. The second lepton in the event must pass the 
same criteria as described above, except for the upper limit on the calorimeter energy within 
a cone associated with its track. The study is performed requiring the reconstructed electrons 
or muons to have P T > 10 GeV. From a total NC sample of 121408 events, 2087 events with 
a second identified electron and 520 events with a reconstructed muon are selected by this 
procedure. Figure la shows the polar angle distribution of the electron with the second highest 
transverse momentum and figure lb shows the polar angle distribution of reconstructed muons. 
The distributions are described by the simulation within the uncertainties, demonstrating that 
the misidentification of a particle as an electron or muon is well understood. 
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The hadronic final state (HFS) is measured by combining calorimeter energy deposits with 
low momentum tracks as described in [19]. Identified isolated electrons or muons are excluded 
from the HFS. The calibration of the hadronic energy scale is made by comparing the transverse 
momentum of the precisely measured scattered electron to that of the HFS in a large NC event 
sample. The transverse momentum of the hadronic system is: 

• P* , which includes all reconstructed particles apart from identified isolated leptons. 

The isolation of identified leptons with respect to jets or other tracks in the event is quantified 
using: 

• their distance Dj et from the axis of the closest hadronic jet in 77— space. For this purpose 
jets, excluding identified leptons, are reconstructed using an inclusive k T algorithm [22- 
24] and are required to have transverse momentum greater than 5 GeV. If there is no 
such jet in the event, D jet is defined with respect to the polar and azimuthal angles of the 
hadronic final state; 

• their distance D track from the closest track in r]-(f) space, where all tracks with a polar 
angle greater than 10° and transverse momentum greater than 0.15 GeV are considered. 

The following quantities are sensitive to the presence of high energy undetected particles 
and/or can be used to reduce the main background contributions. 

• P|, al °, the net transverse momentum measured from all energy deposits recorded in the 
calorimeter. 

• P™ 1SS , the total missing transverse momentum reconstructed from all observed particles 
(electrons, muons and hadrons). P™ 1SS differs most from P£ al ° in the case of events with 
muons, since they deposit little energy in the calorimeter. 

• -^ E , a measure of the azimuthal balance of the event. It is defined as the ratio of the 

Vp 

anti-parallel to parallel components of the measured calorimetric transverse momentum, 
with respect to the direction of the calorimetric transverse momentum [19]. Events with 
one or more high p T particles that do not deposit much energy in the calorimeter (/i, v) 
generally have low values of ^ • 

• <5miss = 2E e — ^ Ei(l — cos 9i), where and 0i denote the energy and polar angle of 
each particle in the event detected in the main detector (9 e < 176°) and E e is the electron 
beam energy. For an event where only momentum in the proton direction is undetected 
5 miss is zero. 

• A02_x, the difference in azimuthal angle between the lepton and the direction of P*. 
NC events typically have values of A0;_x close to 180°. 

• C e 2 = 4PX cos 2 9J2, where E' e is the energy of the final state electron. For NC events, 
where the scattered electron is identified as the isolated high transverse momentum elec- 
tron, Q is equal to the four momentum transfer squared Q 2 . Since the NC cross section 
falls steeply with Q 2 , these events generally have small values of Q. Conversely, elec- 
trons from W decay generally have high values of ( 2 . 
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5 Selection Criteria 



The published HI observation [2] using 1994-1997 e + p data was based on the selection of a 
sample of events with P£ al ° > 25 GeV. This experimental cut mainly selected charged current 
events in a phase space where the trigger efficiency is high. In the selected events all isolated 
charged tracks with transverse momentum above 10 GeV were identified as electrons or muons. 

In the present paper the P£ al ° cut has been lowered to 12 GeV, taking advantage of the im- 
proved understanding of trigger efficiencies with increased luminosity and more sophisticated 
background rejection. The analysis extends the phase space towards lower missing transverse 
momentum for the electron channel (P£ al ° ~ P™ 1SS ) and towards lower P* for the muon chan- 
nel (P|, al ° ~ Pjf)- The lepton identification has also been improved and extended in the forward 
direction. The increased phase space and increased luminosity allow the comparison with the 
SM predictions to be made differentially and with improved precision. Further details of the 
analysis can be found in [25,26]. 

The selection criteria for both channels are summarised in table 1. The dominant back- 
ground in the electron channel is due to NC and CC events. To reduce the NC background, 
events with NC topology (azimuthally balanced, with the lepton and the hadronic system back- 
to-back in the transverse plane) are rejected. For low values of P£ al °, where the NC background 
is largest, a requirement on Q is imposed. A requirement that the lepton candidate be isolated 
from the hadronic final state is imposed to reject CC events. Events which have, in addition to 
an isolated electron, one or more isolated muons are not considered in the electron channel, but 
may contribute in the muon channel. The dominant backgrounds in the muon channel are in- 
elastic muon pair production and CC or photoproduction events which contain a reconstructed 
muon. The final muon sample is selected by rejecting azimuthally balanced events and events 
where more than one muon is observed. 

Following the selection criteria described above, the overall efficiency to select SM W — > ev 
events is 41% and to select SM W — > events is 14%. The main difference in efficiency 
between the two channels is due to the cut on P£ al °, which for muon events acts as a cut on P* 
because the muon deposits little energy in the calorimeter. There is thus almost no efficiency in 
the muon channel for P* < 12 GeV. For values of P* > 25 GeV the efficiencies of the two 
channels are compatible at ~ 40%. 

6 Background Studies 

To verify that the backgrounds (see section 2) that contribute to the two channels are well under- 
stood, alternative event samples, each enriched in one of the important background processes, 
are compared with simulations. For both channels these event samples have the same basic 
phase space definition (9i, Pf , P£ al °) as the main analysis. It should be noted that these selec- 
tions do not explicitly reject signal events, which may be present in the enriched samples. 

The two background enriched samples in the electron channel, defined within the phase 
space 5° < 6 e < 140°, P| > 10 GeV and P| al ° > 12 GeV, are selected with the following 
additional requirements. 
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Electron 


Muon 


Oi 


5° < 9i < 140° 


pi 

T 


> 10 GeV 


pcalo 
r T 


> 12 GeV 


pmiss 

r T 


> 12 GeV 


ryX 
P T 




> 12 GeV 




> 1.0 


^ track 


> 0.5 for 9 e > 45° 


> 0.5 


Q 


> 5000 GeV 2 for P| al ° < 25 GeV 




v p 


< 0.5 (< 0.15 for < 25 GeV) 


< 0.5 (< 0.15 for P£ al ° < 25 GeV) 




< 160° 


< 170° 


# isolated /i 





1 


^miss 


> 5 GeV t 





if only one e candidate is detected, which has the same charge as the beam lepton. 



Table 1: Selection requirements for the electron and muon channels. 
NC enriched sample 

A NC dominated electron sample is selected by requiring Dj et > 1.0. The events in this 
channel mainly contain genuine electron candidates, but with missing transverse momen- 
tum arising from mismeasurement. 

CC enriched sample 

A CC dominated sample is obtained by rejecting events with an isolated muon and apply- 
ing cuts to suppress the NC component. These criteria are £ 2 > 2500 GeV 2 , y 2 - < 0.15, 
^miss > 5 GeV and A0 e _x < 160°. In this sample the missing transverse momentum is 
genuine, but an electron candidate is usually falsely identified. 

The two samples designed to study the backgrounds in the muon channel, defined within 
the same phase space 5° < 9„ < 140°, P£ > 10 GeV and P£ al ° > 12 GeV, are selected with the 
following additional requirements. 

LP enriched sample 

A sample of events predominantly from the two photon process is selected by requiring 
at least one isolated muon and ^ < 0.2 to suppress photoproduction events. 

CC enriched sample A sample dominated by CC events is selected by requiring ^ < 
0.15 and requiring at least one muon candidate that need not be isolated. This selection 
tests fake or real muons observed in events with genuine missing P T . 

The distributions of all quantities used in these selections are well described in both shape 
and normalisation by the SM expectation in regions where there is little contribution from W 
production. This gives us confidence that the backgrounds are described within the uncertainty. 
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Example distributions of the background enriched event samples for the e + p data are shown in 
figure 2 for the electron channel and in figure 3 for the muon channel. Also included in the 
figures are the SM expectations from all processes together and the signal expectation alone. 
Agreement is also obtained between the data and the simulation in all distributions for the e~p 
data sample. 



7 Systematic Uncertainties 

The systematic uncertainties on quantities which influence the SM expectation and the measured 
cross section (see section 8.1) are presented in this section and discussed in more detail in 
[19,26]. The uncertainties on the signal expectation and the acceptance used in the cross section 
calculation are determined by varying experimental quantities by ± 1 standard deviation and 
recalculating the cross section or expectation. The experimental uncertainties are listed below 
and the corresponding variation of the cross section is given in table 2. 

• Leptonic quantities 

The uncertainties on the 9i and the <\>i measurements are 3 mrad and 1 mrad respectively. 
The electron energy scale uncertainty is 3%. The muon energy scale uncertainty is 5%. 

• Hadronic quantities 

The uncertainties on the 9 and <\> measurements of the hadronic final state are both 20 
mrad. The hadronic energy scale uncertainty is 4%. The error on the measurement of X^- 
is±0.02. 

• Triggering / Identification 

The electron finding efficiency has an uncertainty of 2%. The muon finding efficiency has 
an error of 5% in the central (6^ > 12.5°) region and 15% in the forward < 12.5°) 
region. The uncertainty on the track reconstruction efficiency is 3%. The uncertainty on 
the trigger efficiency for the muon channel varies from 16% at P* = 12 GeV to 2% at 
P* > 40 GeV. 

• Luminosity 

The luminosity measurement has an uncertainty of 1.5%. 

• Model 

A 10% uncertainty on the model dependence of the acceptance is estimated by comparing 
the results obtained with two further generators which produce W bosons with different 
kinematic distributions from those of EPVEC. The generators used are an implementation 
of W production within PYTHIA and ANOTOP, an "anomalous top production" gener- 
ator, using the matrix elements of the complete process e + q^e + t^e + b + W as 
obtained from the CompHEP program [27]. 

Contributions from background processes, modelled using RAPGAP, DJANGO and GRAPE, 
are attributed 30% systematic errors determined from the level of agreement observed between 
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the simulations and the control samples (see section 6). The uncertainties associated with lep- 
ton misidentification and the production of fake missing transverse momentum are included in 
these errors. 

A theoretical uncertainty of 15% is quoted for the predicted contributions from signal pro- 
cesses (predominantly SM W production). This is due mainly to uncertainties in the parton 
distribution functions and the scales at which the calculation is performed [5]. 



Source of systematic uncertainty 


Error on measured cross section 


P£ < 25 GeV 


P£ > 25 GeV 


Leptonic quantities 


±0.6% 


±0.6% 


Hadronic quantities 


±3.3% 


±8.3% 


Triggering/Identification 


±3.7% 


±4.7% 


Luminosity 


±1.5% 


±1.5% 


Model Uncertainty 


±10% 


±10% 



Table 2: Summary of experimental systematic errors on the measured cross section in two 
regions of P* . 



8 Results 

For the e~p data sample one event is observed in the electron channel. The kinematics of the 
event are listed in table 3. No events are observed in the muon channel. This compares well to 
the SM expectations of 1.69 ± 0.22 events in the electron channel and 0.37 ± 0.06 in the muon 
channel. 

In the e + p data sample 10 candidate events are observed in the electron channel compared to 
7.2 ±1.2 expected from signal processes and 2.68 ± 0.49 from background sources. One candi- 
date event in the electron channel is observed to contain an e~. This event was first reported and 
discussed in [2]. Four of the other candidate events contain an e + . The charges of the electrons 
in the remaining five events are unmeasured since the electrons are produced at low polar angles 
and they shower in material in the tracking detectors. In the muon channel 8 candidate events 
are observed compared to 2.23 ± 0.43 expected from signal processes and 0.33 ± 0.08 from 
background sources. Four of the muon events observed in the e + p data sample are among those 
first reported and discussed in [2]. The event discussed in [1] is rejected from this analysis by 
the azimuthal difference (A0 M _x) cut. Four of the events have a positively charged muon, three 
have a negative muon and in one event the charge is not determined. 

Distributions of the selected events in lepton polar angle, azimuthal difference, transverse 
mass and P* are shown in figure 4. The lepton-neutrino transverse mass is defined as 

M T = yj (P™ iss + Pj ) 2 - (P™ iss + Pf f, (1) 

where P™ 1SS and Pf are the vectors of the missing transverse momentum and isolated lepton 
respectively. The figure shows the electron and muon channels combined. Also included is the 
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expectation of the Standard Model. The events generally have low values of lepton polar angle 
and are consistent with a flat distribution in azimuthal difference A0;_x, in agreement with the 
expectation. The distribution of the events in M T is compatible with the Jacobian peak expected 
from W production. The kinematics of the events with P^f > 25 GeV are detailed in table 3. 

In three of the eighteen events a further electron is detected in the main detector (9 e < 176°). 
Taking this to be the scattered electron and assuming that there is only one neutrino in the 
final state and there is no initial state QED radiation, the lepton-neutrino mass M Xv can be 
reconstructed. All three events yield masses that are consistent with the W mass, having values 
of 86lg, 73^7 and 79^ GeV. The observation of a second electron in these three events is 
compatible with the expectation from SM W production, where approximately 25% of events 
have a scattered electron in the acceptance range of the main detector. 

Details of the event yields from the e + p data sample as a function of the transverse mo- 
mentum of the hadronic final state, Pjf , are given in table 4 and 5 for the electron and muon 
channels respectively. The combined results for the electron and muon channels are given in 
table 6. At P* < 25 GeV eight events are seen, in agreement with the expectation from the 
Standard Model. At P* > 25 GeV ten events are seen, six of which have P* > 40 GeV, where 
the signal expectation is very low. The probability for the SM expectation to fluctuate to the 
observed number of events or more is 0.10 for the full P^f range, 0.0015 for P* > 25 GeV and 
0.0012 for P* > 40 GeV. The uncertainties on the SM predictions are taken into account in 
calculating these probabilities. 

An excess is observed at P* > 25 GeV in both sets of e + p data. In the 1994-1997 data 
4 events are observed compared to an expectation of 0.80 ± 0.14. In the 1999-2000 data 6 
events are observed compared to an expectation of 2.12 ± 0.36. 

The method published in [2] has been applied to the 1999-2000 data sample. Using this 
method an excess of events is also seen at P* > 25 GeV in this new data sample: 5 events 
are observed for 2.34 ± 0.29 expected. These 5 events selected by the method of the previously 
published analysis are also found by the analysis presented in this paper. 



8.1 Cross Section 

The observed number of events in the e + p data sample is corrected for acceptance and detector 
effects to obtain a cross section for all processes yielding genuine isolated electrons or muons 
and missing transverse momentum. This is defined for the kinematic region 5° < 9i < 140°, 
P l T > 10 GeV, P^ iss > 12 GeV and D jet > 1.0 at a centre of mass energy 3 of ^ = 312 GeV. 
The definition of isolated electrons or muons includes those from leptonic tau decay. The gen- 
erator EPVEC is used to calculate the detector acceptance A for this region of phase space. The 
acceptance accounts for trigger and detection efficiencies and migrations. The cross section is 
thus 

(N data - N bgd ) 



CA 



(2) 



3 Assuming a linear dependence of the cross section on the proton beam energy. 
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where N data is the number of events observed, N bgd is the number of events expected from 
processes treated here as background (see section 2) and C is the integrated luminosity of the 
data sample. 

The cross section integrated over the full range is 

a = 0.31 ±0.10 ±0.04 pb, (3) 

where the first error is statistical and the second is systematic (calculated as described in sec- 
tion 7). 

This result is compatible with the SM signal expectation of 0.237 ± 0.036 pb, dominated 
by the process ep — > eWX, calculated at NLO [5,7]. The small signal components from 
ep — > uWX and Z production are calculated with EPVEC [6] as explained in section 2. The 
cross section is presented in table 7 split into the regions P* < 25 GeV and P* > 25 GeV. 
Whilst the cross section in the low P* region agrees within errors with the prediction, in the 
high P* region it exceeds the expectation. Table 7 also includes two signal calculations in 
which all components are calculated at LO [5,6]. The calculation in [6] is the default calculation 
implemented in the event generator EPVEC. All the calculations agree within the uncertainties. 



9 Search for W Production in the Hadronic Decay Channel 

Since the dominant SM process that produces events with isolated charged leptons and miss- 
ing transverse momentum is W production, it is interesting to search for W bosons decaying 
hadronically. The search for hadronic W decays is performed using events with two high trans- 
verse momentum jets in 117.3 pb -1 of e + p and e~p data from the period 1995-2000. 

Events are selected with at least two hadronic jets, reconstructed using an inclusive kx 
algorithm, with a transverse momentum P T greater than 25 GeV for the leading jet and greater 
than 20 GeV for the second highest P T jet. The minimum P T of any further jet considered in the 
event is set to 5 GeV. The pseudorapidity rj of each jet is restricted to the range —0.5 < rj < 2.5. 
The dijet combination with invariant mass Mjj closest to the W mass is selected as the W 
candidate. The resolution of the reconstructed W mass is approximately 5 GeV. is defined 
as the transverse momentum of the hadronic system after excluding the W candidate jets. 

A cut on the missing transverse momentum P™ 1SS < 20 GeV is applied to reject CC events 
and non-ep scattering background. NC events where the electron is misidentified as a jet are 
rejected [28,29]. The final selection is made with the cuts Mjj > 70 GeV and |cos6>| < 0.6, 
where 9 is the decay angle in the W rest frame, with the W flight direction in the laboratory 
frame taken as the quantisation axis. This phase space is chosen to optimise the acceptance 
for W events and reduce other SM contributions. The overall selection efficiency for SM W 
production is 43% and is 29% for P* > 40 GeV. 

The main physics background to this search is the production of jets via hard partonic scat- 
tering, which is modelled by PYTHIA and RAPGAP for the photoproduction and deep inelastic 
regimes respectively. The predicted cross section is increased by a factor of 1.2 in order to match 
the observed number of events outside the signal region. 
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The systematic uncertainty on the background prediction includes parton distribution func- 
tion uncertainties, the uncertainty on the jet energy scale and uncertainties due to the misiden- 
tification of an electron as a jet. In quadratic sum these give a total systematic error on the 
background prediction of 23% [28, 30]. The SM W production rate has a theoretical error of 
15%, which is added in quadrature to the experimental uncertainties, resulting in an overall 
error of 21%. 

The Mjj distribution (without the Mjj cut) and the P* distribution (with all cuts) of the 
selected data are compared to the Standard Model in figure 5. The final data selected show 
overall agreement with the SM expectation up to the highest P* values. At P* > 25 GeV, 
126 events are observed compared to 162 ± 36 expected with 5.3 ± 1.1 from W production. 
The expectation is dominated by QCD multi-jet production. For P* > 40 GeV 27 events are 
observed in the data, compatible with the expectation of 30.9 ± 6.7, where the W contribution 
amounts to 1.9 ± 0.4 events. Although there is increasing sensitivity to W production with 
increasing , it is at present not possible to conclude from the hadronic channel whether the 
observed excess of events with an isolated electron or muon with missing transverse momentum 
at high P* is due to W production. 



10 Summary 

A search for events with isolated electrons or muons and missing transverse momentum has 
been performed in e + p and e~p data, using the complete HERA I (1994-2000) data sample. 
The selection has been optimised to increase the acceptance for W production events and it 
extends to lower values of hadronic transverse momentum P* than in previous publications. 

One electron event and no muon events are observed in the e~p data, consistent with the 
expectations of 1.69 ± 0.22 and 0.37 ± 0.06 for the electron and muon channels respectively 
in this relatively low luminosity data sample. In the e + p data sample 10 events are observed in 
the electron channel and 8 in the muon channel. These events are kinematically consistent with 
W production. The expected numbers of events from the Standard Model are 9.9 ± 1.3 and 
2.55 ± 0.44 for the electron and muon channels respectively. 

At low Pjf, the number of observed events in both channels is consistent with the ex- 
pectation. At P* > 25 GeV, however, the 10 observed events exceed the SM prediction of 
2.92 ± 0.49. An excess of events is observed in both the 1994-1997 and the 1999-2000 e + p 
data samples. The observed events are used to make a measurement of the cross section for 
all processes producing isolated electrons or muons and missing transverse momentum in the 
kinematic region studied. 

In a separate search for hadronic W decays, agreement with the SM expectation is found 
up to the highest P* values. The high background in this channel, however, does not allow 
one to conclude whether the excess of isolated leptons with missing P T at high P* is due to W 
production. 
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Table 3: Kinematics and lepton charges of the events at high (> 25 GeV). The invariant 
mass M\ v is only calculated for those events with an observed scattered electron. The signifi- 
cance of the charge measurement in numbers of standard deviations is given in brackets after 
the sign. The first event listed was observed in e~p data. The rest were observed in e + p data. 



Electron 


HI Data 


SM expectation 


SM Signal 


Other SM processes 


Pfi < 12 GeV 


5 


6.40 ± 0.79 


4.45 ± 0.70 


1.95 ±0.36 


12 < P* < 25 GeV 


1 


1.96 ±0.27 


1.45 ±0.24 


0.51 ±0.12 


25 < P£ < 40 GeV 


1 


0.95 ±0.14 


0.82 ±0.13 


0.13 ±0.04 


Pfi > 40 GeV 


3 


0.54 ±0.11 


0.45 ±0.11 


0.09 ± 0.04 



Table 4: Observed and predicted numbers of events in the electron channel for all e + p data. 



Muon 


HI Data 


SM expectation 


SM Signal 


Other SM processes 


12 < P* < 25 GeV 


2 


1.11 ±0.19 


0.94 ±0.18 


0.17 ±0.05 


25 < Pjf < 40 GeV 


3 


0.89 ±0.14 


0.77 ±0.14 


0.12 ±0.03 


P$ > 40 GeV 


3 


0.55 ±0.12 


0.51 ±0.12 


0.04 ± 0.01 



Table 5: Observed and predicted numbers of events in the muon channel for all e + p data. 
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Electron and Muon 


HI Data 


SM expectation 


SM Signal 


Other SM processes 


P£ < 12 GeV 


5 


6.40 ± 0.79 


4.45 ± 0.70 


1.95 ±0.36 


12 < P£ < 25 GeV 


3 


3.08 ± 0.43 


2.40 ± 0.40 


0.68 ±0.14 


25 < P£ < 40 GeV 


4 


1.83 ±0.27 


1.59 ±0.26 


0.24 ± 0.06 


P£ > 40 GeV 


6 


1.08 ±0.22 


0.96 ± 0.22 


0.12 ±0.04 



Table 6: Observed and predicted numbers of events in the electron and muon channels combined 
for all e + p data. Only the electron channel contributes for Pjf < 12 GeV. 





Cross Section / pb 




Measured 


SMNLO 


SM LO 
Diener et al. 


SM LO 
Baur et al. 


P# < 25 GeV 


0.146 ±0.081 ±0.022 


0.194 ±0.029 


0.147 ±0.044 


0.197 ±0.059 


P# > 25 GeV 


0.164 ±0.054 ±0.023 


0.043 ± 0.007 


0.041 ±0.012 


0.049 ± 0.015 



Table 7: The measured cross section for events with an isolated high energy electron or muon 
with missing transverse momentum. The cross sections are calculated in the kinematic region: 



5° < 9i < 140 c 



pi 

T 



> 



10 GeV; P^ 1SS 



> 12 GeV and Dj et > 1.0. Also shown are the 



signal expectations from the Standard Model where the dominant contribution ep — > eWX is 
calculated at next to leading order (SM NLO) [5,7] and at leading order (SM LO) [5] and [6]. 
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Figure 1 : Distributions of lepton polar angle of a second reconstructed electron (a) and a recon- 
structed muon (b) for NC events. The combined SM expectation is shown as an open histogram. 
The total error on the SM expectation (see section 7) is given by the shaded band. 
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Figure 2: The e + p data selected in the NC enriched sample (a,b) and in the CC enriched sample 
(c,d) in the electron channel compared with the combined SM expectation (open histogram). 
The total error on the SM expectation is given by the shaded band. The "signal" component of 
the SM expectation is given by the hatched histogram. N data is the total number of data events 
observed for each sample. N SM is the total SM expectation. 
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Figure 3: The e + p data selected in the LP enriched sample (a,b) and in the CC enriched sample 
(c,d) in the muon channel compared with the combined SM expectation (open histogram). The 
total error on the SM expectation is given by the shaded band. The "signal" component of the 
SM expectation is given by the hatched histogram. N data is the total number of data events 
observed for each sample. N SM is the total SM expectation. 
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Figure 4: The final e + p data selection in the electron and muon channels combined compared 
with the SM expectation (open histogram). The total error on the SM expectation is given 
by the shaded band. The "signal" component of the SM expectation is given by the hatched 
histogram. N data is the total number of data events observed for each sample. N SM is the total 
SM expectation. 
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Figure 5: The dijet mass distribution Mjj (a) and the distribution for Mjj > 70 GeV 
(b) compared with the SM expectation (open histogram) in the W hadronic decay channel 
search. The total error on the SM expectation is given by the shaded band. The W production 
component of the SM expectation is given by the hatched histogram. N data is the total number 
of data events observed for each sample. N S m is the total SM expectation. 
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